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Abstract: In this note we stndy the effects of a magnetic field on transport nsing 
holographic models with broken translational invariance. We show that, after carefnlly 
snbtracting off non-trivial magnetisation cnrrents, it is possible to express the DC 
transport cnrrents of the bonndary theory in terms of properties of a black hole horizon. 
This allows ns to obtain simple analytic expressions for the electrical, thermoelectric 
and heat condnctivity tensors. Onr resnlts apply to both isotropic and anisotropic 
models, inclnding holographic Q-lattices and to certain theories where translational 
invariance is broken by linear sonrces for axions. 
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1. Introduction 

Despite the obvious relevance to understanding the rich phenomenology of strange 
metals, it remains an extremely challenging task to calculate the thermoelectric re¬ 
sponse of a strongly interacting theory. One avenue to make progress is to study 
quantum critical theories, where the linear response coefficients are constrained to take 
a scaling form as a function of temperature, T [1, 2, 3, 4], 

However for theories with a net charge density, p, it is much harder to make progress. 
Indeed, in order to obtain hnite transport coefficients it becomes necessary to introduce 
some mechanism for dissipating momentum. Rather than being intrinsic properties of a 
critical theory, the transport coefficients now depend on the details of how translational 
invariance, for instance, is broken. 

Nevertheless, it is still possible to make progress by working perturbatively in some 
small parameter. In particular, detailed results for transport coefficients have been 
derived both within a model of relativistic hydrodynamics (perturbative in p/T^) [5] 
and using the memory matrix formalism (perturbative in the strength of momentum 
dissipation) [6, 7]. 
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An alternative approach is provided by holography. The last couple of years has 
seen a large amount of progress in obtaining analytic expressions for DC transport 
in holographic models in which momentum conservation is violated in some manner. 
These techniques, which originated in the study of the electrical conductivity in massive 
gravity [8, 9, 10], have subsequently been generalised to lattice models [11, 12, 13], 
theories in which translational invariance is broken by linear axions [14], and to the 
calculation of thermoelectric [15, 16, 17] and Hall conductivities [18]. 

The key advantage of these techniques over other approaches is that, rather than 
being valid only in some perturbative regime, it is possible to obtain exact expressions 
for DC transport. The results are therefore valid at all temperatures (that is, including 
-C p) and for any strength of momentum dissipation. In particular, obtaining these 
expressions for electric Hall transport recently allowed a novel mechanism for obtaining 
an anomalous scaling of the Hall angle to be identihed [18]. 

In addition to this behaviour of the Hall angle, many further anomalous aspects of 
the strange metal transport are evident in the effects of a magnetic held on charge and 
heat transport. The original hydrodynamic approach of [5] to magnetotransport was 
motivated by the large Nernst signals detected by Ong et al [19]. Similarly, unusual 
scaling laws are found in the thermopower, magnetoresistance and Hall Lorentz ratio 
[4]. Motivated by these results, in this note we calculate the full set of DC magne¬ 
tothermoelectric transport coefficients for a large class of holographic models. 

As is typical in holography, the key to performing these calculations is to identify 
radially independent quantities in the bulk that can be identihed with the boundary 
currents [20, 9, 13]. However, we will see that the existence of non-trivial magnetisa¬ 
tion currents complicates the usual discussion [21]. In order for us to obtain radially 
independent quantities, it will be necessary to hrst subtract oh the contribution of 
the magnetisation current. Nevertheless, the end result is that we will still be able to 
express the DC transport currents, and hence response coefficients, solely in terms of 
properties of a black hole horizon. 

In Section 2 we present the details of our holographic models and the calculation of 
the DC magnetothermoelectric transport. In Section 3 we close with a brief discussion 
of the signihcance of these results in the wider context of other approaches to mag¬ 
netotransport [5, 22]. In order to improve the readability of our discussion we have 
relegated certain technical details, such as the dehnition of the energy magnetisation 
density and the results for anisotropic theories, to several appendices. 
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We note that whilst this manuscript was in preparation a paper calculating the mag¬ 
netothermoelectric response in holographic models with massive gravitons appeared 
[22], Whilst the lattice and axion models we are considering here are more general, 
the close connection between these models and massive gravity [11, 14] means that our 
calculations and results for DC transport take a similar form to those in [22], 

2. Thermoelectric Transport in a Magnetic Field 

In this section our goal is to calculate the transport coefficients of simple holographic 
models in the presence of a magnetic held. In particular, we wish to obtain the ther¬ 
moelectric linear response of our theories in response to an applied electric held, E, 
and thermal gradient, VT. 

As has been discussed at length in [5, 21] there are subtleties with dehning these 
quantities in the presence of a quantising magnetic held. In particular, the electric 
and heat currents receive additional contributions from spatial variations in the local 
magnetisation. These additional magnetisation currents must be subtracted out of the 
total current in order to obtain the physical transport currents that couple to external 
probes. We therefore need to decompose the total electric, and heat, (5*-*°*^ 

currents as 


j(tot) = J _|_ j(mag) 

= Q + Q(mag) (2.1) 


The magnetisation currents, and have been studied in detail in [21]. There 

it was shown that, at the level of linear response, the magnetisation currents induced 
by electric and thermal gradients are given by 


hmag) 


Q<“‘> 




T 


( 2 . 2 ) 

(2.3) 


where eij is the 2-dimensional antisymmetric tensor (with ey^ = 1) and fi is the chemical 
potential. The other quantities appearing in (2.3) are the magnetisation density M 
and the energy magnetisation density Me- We will not need a precise dehnition of the 
energy magnetisation density in the main text and so relegate the details to Appendix 

A. 
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The main focus of our attention is the linear response of the physical transport 
currents J and Q. This dehnes the electric, d, electrothermal, d, and heat, k \ con¬ 
ductivities according to 

J = aE — aVT 

Q = aTE - kVT (2.4) 


In the presence of a magnetic field, each of these conductivities are 2 by 2 matrices. 
For isotropic systems, which will be the main focus of our discussion, these can be 
decomposed into their symmetric, e.g. a^x-, and antisymmetric, e.g. axy, parts^. This 
means that the linear response is described by six functions axx-, o'xy, Oixx, ctxy, ^xx, ^xy 


Our goal then, in this section, is to calculate the DC limit of these transport coef¬ 
ficients for a large class of holographic models. Here, we focus on holographic models 
that break translational invariance but preserve the homogeneity of the bulk action. 
These models have received a large amount of recent attention [12, 13, 14, 18, 15], 
following the realisation that is possible to obtain analytic expressions for their DC 
transport properties [9, 13]. For concreteness we will work with the following simple 
bulk action 




d" 


X 




R 


+ ■i>(0)((Sxi)" + (dx-^n + U(0) 




(2,5) 


although our considerations can be trivially extend to more general actions. The trans¬ 
lational invariance of the 2-1-1 dimensional boundary theory is broken by constructing 
background solutions where we have xi = kix and X 2 = k 2 y- Nevertheless, since only 
derivatives of the Xi fields feature in the action, the system remains homogeneous and 
can be studied using ODEs. For simplicity of presentation we will only consider the 
isotropic case ki = k 2 = k in the main text, although we present results for anisotropic 
systems in Appendix B. 


The class of models described by (2.5) includes many of the theories that have been 
studied in the literature. In particular, if we choose to set 0 = 0, <h = const then 
the X fields correspond to massless axions in the bulk. These are dual to marginal 
operators, in the boundary theory, in which translational invariance is broken by 
a linear source xf'^ = kxi. The DC transport properties of these theories have been 
intensely studied and are directly related to those of massive gravity theories [14]. 

^Note that R is not the true thermal conductivity, k, but rather the thermal conductivity in zero 
electric field. These are related hy k = k — Ta.a~^.a. 

^Note that isotropy implies that (lyy = Oxx and that Oxy = —<Jyx- 
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On the other hand it is also possible to obtain an action of the form (2.5) starting 
from the Q-lattice models introduced in [12], These Q-lattice models break translational 
invariance by introducing an oscillatory lattice through two complex scalar helds Ti ~ 
Q^kx ^ \j /2 ~ The canonical action of these charged scalars can be rewritten in 

the form (2.5), which is convenient for discussing transport, by performing the polar 
decomposition^ 'hi(r) = The held 0(r) therefore corresponds to magnitude 

of these lattices in the bulk, whilst the Xi(r) can be thought of as their phase. If one 
uses conventional kinetic terms this results in the action (2.5) with the choice ‘h(0) = 0^ 
and the requirement that we should identify the helds Xi under shifts of Xi —t y* + 27r. 

For isotropic solutions a suitable ansatz for the background metric and gauge held 
takes the form 


ds^ =—Udt^ + U ^dr"^ + {dx^ + dy^) (2.6) 

A = a{r)dt — Bydx (2.7) 

where we assume that the geometry approaches AdS at the boundary r —)■ cxd. The 
temporal gauge held a(r) asymptotes to a constant value y which is interpreted as the 
chemical potential of the boundary theory. Likewise B corresponds to the magnetic 
held in the dual theory. 

In addition we will assume that there is a regular black hole horizon located at 
a position r_|_ in the bulk. Near this horizon we can expand the radially dependent 
background helds as 


U ~ 47rT(r — r+) + ... 

Qj a+(r - r+) + ... 

~ 14 + ... 

(f) ~ 0_i_ + ... (2-8) 

where T is identihed with the temperature of the dual theory. 

As is standard, the transport coefficients are computed in holography by studying 
perturbations of the background solution. We will follow the approach introduced in 

^Note that in the main text we are assuming that the two fields are related by a bulk Z2 symmetry 
and so can have the same radial profile. We allow for anisotropic configurations in Appendix B. 
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[13] and calculate the DC conductivity by applying linear sources to the boundary 
helds. That is we consider the perturbation ansatz 

= -By + {-E + ^a(r))t + SA^{r) 

Ay = 5Ay(r) 

gtx = -itu + e^^5htx{r) 

gty = e^^6hty{r) 

Qrx 6 Shrxij'^ 

Qry ^ 

Xi = kx + 6xiir) 

X2 = ky + 6x2ir) (2.9) 

which corresponds to applying an external electric held Ei = E6ix and temperature 
gradient (VT), = ^SixT to the boundary theory. 

Electrical Currents 

As was hrst realised in [9, 11, 13], the key reason it is possible to calculate transport 
coefficients in these models is that the currents of the boundary theory can be related 
to radially-independent quantities in the bulk. In particular, the AdS/CFT dictionary 
tells us that the expectation value of the currents are given by the quantities 

as r^oo (2.10) 

Usually the huxes ^/^Z{(f))E^'' are independent of radial position, and so can be 
evaluated anywhere in the bulk. Evaluating these huxes at the horizon allows the 
conductivity tensor to be extracted. 

However, in this case, the presence of non-trivial magnetisation currents complicates 
the discussion. It is simple to use the linearised Maxwell equation dy{yJ—gZ{(j))E'^^^) = 0 
to show that for the perturbations (2.9) we have 

dr{^Z{cl>)En = -dt{^/^Z{<p)E^^) = 0 

driV^Z{(j))Ey'') = -d,{^Z{ct>)Ey^) = -e-^^Z{ct>)Bi ( 2 . 11 ) 

which implies that the huxes are no longer constant in the presence of a thermal gra¬ 
dient. Note that if we had also applied a thermal gradient in the y direction, we would 
also have found that ^J—gZ{(())E^^ depended on the radial coordinate, r. 
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Nevertheless, we can still construct quantities that are independent of the radial 
coordinate by integrating equations (2.11). That is we dehne bulk fluxes by 


J%r) = 

jy(r) = ^Z{(j))Fy^ - iM{r) ( 2 . 12 ) 

where M(r) is dehned to be 

nr 

M{r) = - dfe-^^Z{4))B (2.13) 

J r+ 

This extra term has been chosen so that the currents JT”* defined in (2.12) are radially 
constant by construction, 


drT = 0 (2.14) 

However, as we approach the boundary, r —)■ cx) we can no longer identify them with the 
total currents in the boundary theory. Nevertheless, to calculate the response coeffi¬ 
cients it is the transport currents, rather than the total currents, that we are interested 
in. Remarkably, we show in Appendix A that as r —)■ cxd then M(r) corresponds to 
the magnetisation density of the boundary theory. The effect of the additional term 
in (2.12) is therefore simply to subtract off the magnetisation current so that near the 
boundary we have 

(/*) = J^r) as r —)■ cxD (2-15) 

i.e. these constant bulk fluxes JT”* precisely correspond to the transport currents of the 
boundary theory^. 

Having related the transport currents to bulk constants we can proceed to calculate 
the DC transport as normal. Linearising these expressions according to (2.9) gives the 
bulk constants 

= -Z{(t))U5A'^ - Z{(j))e^^a'6ht^ - BZ{^)U6hry 
jy = -Z{(t))U5A'y - Z{(j))e^^a'6hty + BZ{^)U6hrx - ^M{r) (2.16) 

Since these expressions are independent of the radial coordinate, we can choose to 
evaluate them wherever we like. The trick, as always, is to proceed to the horizon 

"‘Note the the reason we have only had to subtract off the magnetisation current from is because 
our ansatz (2.9) only corresponds to applying a thermal gradient in the x direction. 
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where the constraints of horizon regularity imply that 

E 

6Ai = -—A-\n{r - r+) + 0(r - r+) 

47rj 

^Xi = 0{{r-r+f) 

Shu = UShri - ^J^^y ln(r - r+) + 0(r - r+) (2.17) 

Note that for our ansatz we have Ei = Si^E and = Six^ but we have left the regularity 
conditions in their general form. 


Furthermore the dehnition of M(r) implies that it vanishes at the horizon. We 
therefore hnd that that the transport currents can be expressed solely in terms of 
properties of the horizon^. In particular we have that 


JX 

jy 


Z{(j))Ex - e^^Z{(^)a'5htx - Z{(j))BShty 


r+ 


Z{cj))Ey - e^^Z{cj))a'Shty + Z{^)B5htx 


r+ 


( 2 , 18 ) 


All that remains is to determine the values of the graviton fluctuations Shu at the 
horizon. This can be done by examining the t — x and t — y components of the linearised 
Einstein equations which read 


U{e*^Sh[J - {B^Z + e^^k‘^^)Shtx + BZUe^^a'SKy = -e^^ZUa'Sa'^ 

Uie^^Sh'tyY - {B^Z + e^^k^^)Shty - BZUe^^a'SKx = -e^^ZUa'Sa'y + BZ{-E + ^a(r)) 

(2.19) 


After imposing the regularity conditions (2.17) these reduce to requiring that we satisfy 

{B^Z{<t)) + e^^k^^{<t)))Shtx - BZ{(j))e^^a'Shty = Z{(j))a'E + 

{B^Z{^) + e^^k‘^^(j)))Shty + BZ{(j))e^^a'Shtx = BZ{(j))E (2.20) 


at the horizon. It is then straightforward to invert these equations and substitute for 
Shu into (2.18). The resulting expressions for the electrical currents can be compared 
to (2.4), which allows us to read off the electrical conductivity tensor as 


(T- 


^xy 


e 2 ^A; 2 $(p 2 + 52^2 ^ ^e2V^2^) 


B2p2 + ( 52 ^ + g2V^2$)2 

(p2 + 52^2 + 2Ze2^A;2<I>) 
^^fiV + (52Z + e2'"A;2<l>)2 


r+ 


r+ 


( 2 . 21 ) 


^This should be contrasted with the total current which, since it depends on the magnetisation, is 
sensitive to the full geometry. 



whilst the electrothermal conductivities are 


= 


spe 


2Vp2^ 


a 


xy 


= sB 




B2p2 + [B‘^z + e‘^^k‘^^f\r+ 


which we have expressed in terms of the boundary charge density p 
entropy density s = 47re^'^|r+. 


( 2 . 22 ) 
—Ze^'^ a' and 


Heat currents 

Up until now we have focused solely on the electrical current. In order to extract the 
heat conductivity, we need to consider the heat currents of the boundary theory. In 
[15] it was shown that this can be done by considering the bulk two-form dehned 
by 

+ ^(2a(r) + Ex)Z{^)F^^’^ (2.23) 

where is the vector held dt- The heat currents can then be identihed with this two 
form in a similar way to how the electrical current is related to the held strength. In 
particular, at the linearised level we can make the identihcation 

as r^oo (2.24) 

which follows from evaluating G” for the perturbations (2.9) to get 

(Qh°th) = G^( ^ - a(r)^^Z(0)F^^ as r ^ oo (2.25) 

Up to contact terms, the hrst term in this expression corresponds to the expectation 
value of the energy momentum tensor (X 5 ]_ xhe second term subtracts oh the 

electric current to get the heat current (5*-*°*^* = T*^*°*)*° — 

Once we have these huxes, much of our earlier discussion can now be applied to 
the heat currents. The motivation for introducing the two-form is that, in the 
absence of a thermal gradient, it was shown in [15] to satisfy d^{yJ—gG^^) = 0. As 
a result the linearised huxes \/—gG^^ were independent of the bulk radial coordinate. 
However, the existence of magnetisation currents means that this is no longer true for 
the perturbations (2.9). Rather we have that 

driV^G^ = -dtiV^G^n - dyiV^Gy^ 

= 0 

dr{y/^G^y) = -dti^G^y) - d^i^G^y) + e-^^Z{(j))B^a{r) 

= -e-^^Z{<t>)B{E-2ia{r)) 
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It is therefore again necessary to add an extra term to the fluxes in order to obtain 
radially independent constants. That is we construct 

Qv = t/2 (—- a{r)^Z{<j))Fy^ - M{r)E - 2MQ{r)i (2.26) 
where Mqlr) is given by 

nr 

Mqir) = / dfe-^^ Z{(j))Ba{f) (2.27) 

J r+ 

The additional terms in the dehnition of ensure that these modihed fluxes are 
radially constant 


= 0 (2.28) 

by construction. Once again the fact we have been forced to introduce this extra term 
reflects the presence of magnetisation currents in the boundary. In Appendix A, we 
show that as r —)■ cxD then Mqlr) precisely approaches the heat magnetisation density, 
Mq = Me — pM, of the dual theory. The effect of this additional term is therefore to 
subtract off the contribution of the magnetisation current from (2.25). More precisely, 
we have that near the boundary 


(Q*) = Q as r^oo (2.29) 

i.e. the bulk constants Q* correspond to the heat transport currents of the boundary 
theory. 

We can now repeat the trick we used with the electrical currents and evaluate these 
constants at the horizon. The dehnitions of Mq{r) and M(r) imply that they vanish 
at r = r_|_ and so, as for the electrical case, we see that the transport currents can be 
expressed locally in terms of horizon helds. Using the regularity conditions (2.17) we 
hnd that these expressions take the simple form 

Q" = 

QV = -u'e^^6hty\r+ (2.30) 
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Fortunately, in our discussion of the electrical current we have already determined 
the values of 6hti\r+ by inverting (2.20). We can therefore extract the thermoelectric 
conductivity, a, as 


(Utt = 


spe^^k'^^ 


a 


xy 


= sB 


B‘^p^ + {B^Z + e‘^^k‘^^)Ar+ 
{p^ + B‘^Z^ + 


B2p2 + + e^^k^<^y\r+ 


(2.31) 


which reassuringly agrees with the expression we obtained from the electrical current 
(2.22). Finally, the heat conductivity, R reads 


^xx 


^xy 


s^T{B‘^Z + e2^A;2d>) 

B2p2 ^ (52^ + g2V^2$)2 

s^TpB 

B2p2 + (B2^ + e2V'A;2$)2 


r+ 


r+ 


(2.32) 


3. Discussion 

Whilst on first glance the expressions for these transport coefficients seem rather 
baroque, on closer inspection they display a remarkable simplicity. Perhaps the most 
striking aspect of the equations is that, just as in the hydrodynamic analysis of [5], the 
entire set of DC transport coefficients are described by two parameters. That is, aside 
from thermodynamic factors, it is only the functions® 

accs = Z{(j))\r+ ^ = e^^k‘^^r+ (3.1) 

T 

that appear in the transport coefficients. Here we have dehned the timescale r so that 
with the above identihcations the electric conductivity tensor takes precisely the same 
form as in [5], where corresponded to the momentum dissipation rate. 

Although the agreement of the electrical conductivity tensor between these two ap¬ 
proaches is striking, it does not extend to the thermoelectric response coefficients (as 
has previously been emphasised in [16, 17, 22]). This is not necessarily a surprise since, 
as we stressed in the introduction, as soon as one includes a net charge density then the 
transport coefficients will depend on the microscopic way in which momentum dissipa¬ 
tion is incorporated. The holographic results we have obtained therefore suggest that 
the mechanism for momentum dissipation in holography is different from the particular 
model studied in [5]. 

®Here £ is the energy density, V the pressure and we have identified Uccs as a ‘charge-conjugation 
symmetric’ conductivity. 
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Whilst our results differ from those presented in [5], they take the same qualitative 
form as the recent results using massive gravity [22], In particular, we reproduce 
the results of linear axions (and hence massive gravity [14]) by choosing to take the 
fnnctions ^(0) and <h(0) to be constant. Althongh the strnctnre of the eqnations is 
very similar, we emphasise that onr results apply for much more general holographic 
models - regardless of the form of these fnnctions. 

In particular, by varying the choice of action, it is now possible to obtain qnite gen¬ 
eral scalings (with temperatnre) in the horizon qnantities Z{(f))\r_^_ and <l>(0)|,._^. It is 
therefore tempting to see if one can choose these scalings to match the phenomenol¬ 
ogy of the cnprates. This idea was anticipated in [22] where, after matching r and 
aces to the Hall angle as proposed in [18], the scalings of the thermoelectric and heat 
transport coefficients, to leading order in B, conld be dedneed. However, since the 
holographic results differ in general from the hydrodynamic analysis, it is not yet clear 
to what extent these scalings are nniversal and hence can be meaningfnlly compared 
to experiment. 

Nevertheless, the models and results we have presented here should be capable of 
realising these proposed scalings explicitly. Furthermore, since onr results are valid 
even in strong magnetic helds, it is possible to use these models to go beyond leading 
order in H. In particular, it has recently been proposed that the magnetoresistance 
of strange metals can be attributed to the effects of the magnetic field on the critical 
theory itself [23]. It would therefore be extremely interesting to use the holographic 
models discussed here to study transport in the presence of a strong magnetic held. 
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A. Magnetisation and Energy Magnetisation Densities 

In this Appendix we wish to derive the formnlae for the magnetisation and energy 
magnetisation densities used in the main text. The dehnition of the magnetisation is 
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familiar. If we apply a magnetic field to the boundary theory via a source = —By, 
then the magnetisation density is given by differentiating the Euclidean action Se, as 


M 


1 OSe 


(A.l) 


where V is the volume of the boundary held theory. 


The energy magnetisation density is dehned as an analogous quantity for the metric. 
That is we should apply a source = —Biy and differentiate with respect to Bi 


Me = 


IdSE 

VdBi 


Bi=0 


(A.2) 


We now wish to calculate these for the background solutions to our action (2.5). To do 
this, it is convenient to consider solutions obeying the ansatz 


X* = kxi 0 = 0(r) 


At = a(r) A^ = -By + (a(r) - y)Biy 


ds^ = —U{r){dt + BiydxY + . + e^^^'^Hdx'^ + dy"^) 

U[r) 

When Bi = 0, these are simply the background solutions we studied in the main text. 
At leading order in Bi we also have applied a source = —Biy to the boundary the¬ 
ory which will allow us to evaluate the energy magnetisation density. The higher order 
terms in Bi in the metric ensure that we have written down a consistent ansatz. In par¬ 
ticular note that even though we have introduced dependence on the ^/-coordinate, this 
does not appear in the equations of motion. The dynamical helds 0(r), a(r), U{r), V{r) 
only depend on the radial coordinate. 


In order to calculate the magnetisation and energy magnetisation densities, we need 
to differentiate the action with respect to B and Bi, before setting Bi to zero. To do 
this we hrst differentiate the off-shell bulk action, before evaluating these derivatives 
on the equations of motion. With our ansatz we hnd that the Einstein-Hilbert term in 
the action (2.5) can be written as 




dr 


e 


-2V 


2 



2e^^(l7" + 417'V' + 617V'^ + 417V") 
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Similarly the scalar terms take the form 


^ cscalar 


dr 


rl 


'r+ 




Finally we have the Maxwell term 
1 


^ cMaxwell 


poo -2V P 

I dr-Z(0) (5 + 5i/i)2-25i(5 + 5i/i)a + 5iV-e^V 

r+ 2 L 


Although we will not need them here, the equations of motion can of course be deduced 
by varying this action with respect to the dynamical helds a(r), 0(r), U{r), V{r). 


To get the magnetisation density, we simply need to set -Bi = 0 and then differentiate 
the action. We note that the only explicit B dependence is in the Maxwell term and 
hence we hnd^ 


M =Are-^'-Z(4,)B (A.3) 

It is now clear to see that the M(r) dehned in (2.13) is precisely the magnetisation 
density when r —)■ cxd. 

Likewise we can construct the energy magnetisation density by first differentiating 
with respect to Bi and then setting i?i = 0. At linear order in Bi we again hnd that it 
is only the Maxwell term that contributes®. We thus read off the energy magnetisation 
density as 


Mr = 


1 dS^ 


VdBi 


dre Z{(j))B{^ — a{r)) 


Bi=0 


(A.4) 


/ r+ 


Finally we deduce that the heat magnetisation density is given by 


Mq = Me- fiM = I Z{(j))Ba{r) 


(A.5) 


' r“+ 


from which we can see that the function Mqlr) (2.27) is equivalent to Mq when r —)■ cxd. 


^Note that the Euclidean action constructed via a Wick rotation t —>■ —ir has an extra minus sign 
relative to the Lorentzian action (2.5). 

®Note that at leading order in Bi the magnetisation and energy magnetisation densities do not 
receive any contribution from the boundary counterterms. 
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B. Anisotropic Models 

It is straightforward to generalise our calculations to anisotropic theories. In particular, 
we can consider the action 

s = J R - ^[{d(j)Y + + $ 2 ( 0 ) (<9X2)^] + l^r( 0 ) - 

where we now break translational invariance in the x and y directions by constructing 
background solutions with xi = kix and X 2 = ^ 2 l/- In order to allow for anisotropic 
solutions we should also modify our metric ansatz to 

ds^ = —Udt^ + U~^dr‘^ + e^'^^dx^ + e^^'^dy'^ (B-1) 

The resulting expressions for DC transport are more complicated that the isotropic 
case, but are again simplihed somewhat by introducing the thermodynamic factors 
s = and p = They read 


/T 

eyi+V2^2^2(p2 + 52^2 + Ze2'^2^2^i) 


(B.2) 

^ XX 

“ B^p^ + {B'^Z + e‘^^^kl^i){B‘^Z + e^^^kl^ 2 ) 

r+ 

rr — 

(p2 + 52^2 + Ze2^2A;2$i + Ze2^iA;2$2) 

■ + e‘^^^kld>i){B^Z + 


(B.3) 

^ xy 

) 

r\i 

spe^^+^^kld >2 


(B.4) 

^XX 

~ B^p^ + {B^Z + e‘^^^kl^i){B^Z + 

»■+ 

r\ 

sB{p^ + B^Z^ + Ze2^2^2^i) 


(B.5) 

ULxy 

“ B2p2 + (^2^ + e‘^^^kld>i){B^Z + e2^iA;2$2) 

r+ 

K 

47 re^^^sT{B^Z + e^^^kl^ 2 ) 


(B.6) 

'^XX 

~ B^p^ + {B^Z + e‘^^^kld>i){B^Z + e^^^kl^ 2 ) 

r+ 


s'^TpB 

~ B^p'^ + {B^Z + r-+ 


(B,7) 


where the expressions for the remaining transport coefficients (e.g. (jyy) can be trivially 
obtained from those we have presented above through swapping around the labels 1 
and 2. 
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